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The leading correction to scaling associated with departures of the initial condition from the 
scaling morphology is determined for some soluble models of phase-ordering kinetics. The result 
for the pair correlation function has the form C{r,t) = fo{r/L) + f\{r / L) + • ■ •, where L is a 
characteristic length scale extracted from the energy. The correction-to-scaling exponent ui has the 
value uj = A for the d = 1 Glauber model, the n-vector model with n = oo, and the approximate 
theory of Ohta, Jasnow and Kawasaki. For the approximate Mazenko theory, however, uj has a 
non-trivial value: uj = 3.8836 ■ ■ ■ for d = 2, and u = 3.9030 ■ • • for d = 3. The correction-to-scaling 
functions fi(x) are also calculated. 
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I. INTRODUCTION 

The phenomenon of phase-ordering kinetics has at- 
tracted considerable attention in recent years It 
is now well established that, except in certain excep- 
tional circumstances, a scaling regime develops at late 
times, in which the order-parameter morphology is time- 
independent up to one overall time-dependent length 
scale L(t). This means, for example, that the equal- 
time pair correlation function of the order parameter, 
C{r,t) = ((/)(x -I- r, i)0(x, i)) has the asymptotic scal- 
ing form C{r,t) — f[r/L{t)]. Relatively little attention, 
however, has been devoted to the study of how the scal- 
ing regime is approached, i.e. the form of 'corrections 
to scaling' in phase-ordering kinetics. In particular, an 
understanding of corrections to scaling is important in 
interpreting experimental or simulation data, and in ex- 
tracting asymptotic scaling exponents and functions. 

This article attempts a first systematic study of cor- 
rections to scaling. We show that there is a correction- 
to-scaling exponent lu associated with the deviations of 
the order parameter morphology from the scaling mor- 
phology, i.e. with the fact that the initial state is not in 
general the scaling state. Presupposing a suitable defini- 
tion of L{t), which will require some discussion, the lead- 
ing corrections to scaling in the pair correlation function 
enter in the form C{r,t) — fo{r / L) + L^"^ fi{r /L), where 
/o(x) is the 'scaling function' and fi{x) the 'correction- 
to-scaling function'. There are indications that lo is, in 
general, a nontrivial exponent of phase-ordering kinetics. 
Although for most of the simple, exactly soluble mod- 
els presented here it has the value w = 4, we find that 
it takes nontrivial values within an approximate calcula- 
tion for more realistic models. The models considered are 
the one-dimensional (ID) Glauber model, the 0{n) non- 
linear sigma model for n — oo, the ID time-dependent 
Ginzburg-Landau (TDGL) equation, the Ohta-Jasnow- 
Kawasaki (OJK) approximation ||] for the general TDGL 
equation, and lastly the Mazenko approximation It 
is the last of these that yields non-trivial values for lo. In 
fact LO enters the theory in a similar way to the exponent 
A, which describes the decay of the autocorrelation func- 
tion and is known to be nontrivial B. For the ID TDGL 



equation, we find that the approach to the scaling state is 
much faster than for the other models: for typical initial 
conditions the corrections to scaling vanish exponentially 
fast as a function of the smallest domain size. 

It is important to recognize that corrections to scaling 
can have more than one origin. For a scalar order param- 
eter, for example, in which the morphology is specified 
by the locations of the domain walls, the scaling regime 
requires that the domain scale L(t) be much larger than 
the width ^ of the walls. Therefore one expects a cor- 
rection to scaling associated simply with the non-zero 
width of domain walls, and entering as a power of £,/L. 
This is a distinct contribution from the one considered 
here, which is associated with departures from the scaling 
morphology in the initial condition. This latter contri- 
bution will survive even in the thin wall limit described 
by, for example, the AUen-Cahn theory [Q]^ in which 
infinitely thin interfaces are driven by their local curva- 
ture. Throughout this paper we restrict ourselves to this 
thin- wall limit, or the corresponding 'hard-spin' limit (or 
non-linear sigma model) for a vector order-parameter. 
This leads to useful simplifications. 

The paper is organized as follows. Section II deals with 
the ID Glauber model, and introduces some general con- 
cepts. The large-n vector model, ID TDGL equation, 
OJK Theory, and Mazenko theory, are covered in sec- 
tions III to VI respectively. Section VII concludes with 
a discussion and summary. 



II. THE ONE-DIMENSIONAL GLAUBER MODEL 



In the continuum limit, the equation of motion for the 
pair correlation function C (r, t) has the form |^ 



dtC = dlC, C(0,t) = l. 



(1) 



The constraint C(0, i) = 1 can be built in through a 
source term at r = 0, i.e. we modify the equation of mo- 
tion to 



dtC = dlC + A{t)5{r), 



(2) 
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where A{t) is chosen to satisfy the constraint. Let us 
define L to be the average domain size. Then for r L, 
we have 



C(r,t) = l-2p^|r| = l-2|r|/L, 



(3) 



where pw — ^/ L is the density of domain walls, and we 
have allowed for the possibility of either one or zero walls 
in an interval of length r. Equation (||) is then correct to 
order r. 

Integrating across an infinitesimal interval around 
r = 0, and comparing with (|), gives A — A/L. Putting 
this into (H), and writing C{r,t) in terms of r and L in 
the form 



C{r,t)^ f{r/L,L) 
gives the following equation for f{x,L): 



(4) 



(5) 



where dots and primes indicate derivatives with respect 
to t and X respectively. 

In the limit L — > oo, we expect f{x, L) to approach the 
i-indepcndcnt scaling function fo{x). Balancing powers 
of L on the right of (H) then implies LL — constant, i.e. 
L oc t^/^. Including the leading correction to scaling, of 
relative order L~", we write 



f{x,L) = fo{x) + L-'^h{x) + ^ 



(6) 
(7) 



in (|^), and equate terms of leading {L ^) and next-to- 
leading (i^^^'^) order. This gives 



-I- axf^ + U{x) = 
/f + axfi + atuh + bxf^ = 0. 



(8) 
(9) 



Integrating (||) with the boundary conditions /(O) 1, 
f'{0+) = -2 [from (|)] gives (for a; > 0) 



a = 27r 



/o(x) =erfc( ^^|-x 



(10) 
(11) 



Equation (g) for /i becomes 

f[' + axf[ + aujfi = hx exp(-aa;V2) , (12) 
2b. Writing Ji{x) = 



where b = 6(2a/7r)^/^ 
exp(— 0x^/2)3(2;) gives 



— axg' + a(cj — l)g ~ bx . 



(13) 



What are the boundary conditions on g{x)7 Clearly 
5(0) = 0, because C{0,t) = 1 is already implemented by 
/o(0) = 1. Similarly, the condition 9rC(r, i) I r=o = — 2/i, 
from (0), is already guaranteed by /o(0) = —2, so g'{0) = 



0. Furthermore, ( |T^ ) then gives g"(0) = 0, implying that 
the series expansion for g{x) starts at 0{x^). Inserting 
the series solution g{x) = Y^^=3, 9nx"' gives 53 — 6/6, and 
the recurrence relation gn+2 = [a{n + 1 ~ to) / {n + l){n + 
2)]gn for the higher-order odd coefficients, all even coeffi- 
cients vanishing. In order that /i (x) decrease faster than 
a power-law for large- a;, as required on physical grounds 
for initial conditions with only short-range correlations, 
the series expansion for g{x) must terminate. This gives 
the condition uj — n + 1 — 4,6, 8,---. We conclude 
that the leading correction-to-scaling exponent for the 
1-d Glauber model is 

w = 4 , (14) 
with corresponding correction-to-scaling function 



/i(x) = (6/6) a;^ exp(-aa;V2) • 



(15) 



III. THE LARGE-N LIMIT OF THE 0(N) MODEL 

It is convenient to work with spins of fixed length, 
i.e. with the nonlinear sigma model, to avoid additional 
(uninteresting) corrections to scaling associated with the 
gradual saturation in the length of the spins as the coars- 
ening proceeds. The non-conserved dynamics of the 0( 
non-linear sigma model is described by the equation 



dtcj)^\/^cf>+{S/cf>f 0, 



(16) 



corresponding to the equation dt4> = —SF/Scj), with free 
energy F = d'^ x {\/ (j))'^ , subject to the constraint 

[^(X,t)]2 = l. 

In the limit n ^ 00 we can replace (V0)^ by its mean 
in the usual way 0|. Let us call this mean a{t). In the 
scaling regime, dimensional analysis gives a{t) ~ \/2Lp' , 
where A is a constant. The energy density is just e{t) = 
a{t)/2 = X/AL^. It is, in fact, convenient to define L{t) 
through this relation for all times t (in the same spirit 
as the Id Glauber model). Then multiplying ( p^ by 
(/)(x -I- r, t) and averaging over initial conditions gives the 
equation 



1 dC_ 

2 'dt 



2i2 



C 



(17) 



for the pair correlation function C (r, t) . 

In analogy to our treatment of the Id Glauber model 
we write C{r,t) = f{r/L,L). Then (|l7| ) becomes the 
following equation for f{x, L): 



2 dL V x ■' 2-' 



^^/'. (18) 



Inserting the forms (^ and (^ (taking a = 1/2 without 
loss of generality) , and equating coefhcients of the terms 
in L^^ and as before gives 
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Consider first /q. The two linearly independent solu- 
tions have large- a; behavior ~ x~^^ and ~ exp(— a;^/8). 
The power-law term must be absent for a physical solu- 
tion corresponding to an initial condition with only short- 
range correlations. Setting fo{x) — exp(— (70(2;) in 
® gives <?(,' + [(d-l)/x-x/4]5^,+ (A/2-d/4)go - 0. The 
boundary conditions are 30(0) — Ij 9oiQ) = (the lat- 
ter being required by the differential equation, to avoid 
a singularity at x = 0). To retain the gaussian tail in 
/o, the series solution for go must terminate. This fixes 
A = d/2 -\- n (n = 0, 1, 2, • ■ •), with a corresponding set 
of polynomial solutions for 170 1 the first two of which are 
gj^^ = 1, and ^g^-* = 1 — /Ad. An explicit solution of 
C for a general initial condition shows that these differ- 
ent scaling solutions are selected by the small-Zc behavior 
of the structure factor, iS'(k, t) [the Fourier transform of 

C(r, i)], at t = 0. The polynomial solution g^^^ corre- 
sponds to an initial condition with 5(k, 0) oc (fc^)" for 
fc — > 0. A generic initial condition, therefore, selects the 
n — solution, i.e. A — d/2 and /o(a;) — exp(— 

Consider now the correction-to-scaling function /i. 
Putting fi{x) = exp{-x^/8)gi{x) in @, with A = d/2 
and /o(a;) — exp(— gives 

Again one seeks a series solution for gi{x). The bound- 
ary conditions .gi(0) = = 5'i(0) imply that the sc- 
ries has the form .91(2;) = J2^=2 9nx"- Substituting 
this into the differential equation one readily finds that 
gn = for n odd, g2 = 0, 34 = b/[32{d + 2)], and 
= [{n — t^) /4(n -I- 2) (n -f d)] <?„ for even n > 4. For a 
physically sensible large-x behavior, the series must ter- 
minate, just as in the one-dimensional Glauber model. 
This requires oj — 4,6,8,---. For each such lu, there 
is a corresponding correction-to-scaling function. The 
generic case, corresponding to short-range correlations 
in the initial condition, is the smallest value of lu, i.e. 
LU = 4, with gi{x) = bx'^ /32(d + 2). The correction-to- 
scaling function is then 

^^(^)-32(dT^^^H-^) • ^''^ 

The value of uu for the n = oa limit, cj = 4, is thus 
identical to that of the \D Glauber model, raising the 
question of whether this could be a general result for 
models with non-conserved dynamics. Indeed, the Ohta- 
Jasnow-Kawasaki theory discussed in section V gives the 
same result. The related approach of Mazenko (section 
VI), however, gives a different, and indeed nontrivial, re- 
sult. In the following section, we discuss another model 
with a; 7^ 4. 



IV. THE ONE-DIMENSIONAL TDGL EQUATION 

Another exactly soluble model (in a suitable limit) is 
the time-dependent Ginzburg-Landau (or TDGL) equa- 
tion in one dimension. This reads 

dt^ = dl<t> - V'(0), (23) 

where 4>{x, t) is the order-parameter field and V^(</>) is a 
potential function, with a symmetric double-well struc- 
ture. For convenience, we can take the minima of V to 
be at = ±1. Eq. ( [2^ ) represents the simple relaxational 
dynamics dt4> = —SF/S(j), with free-energy functional 
= J dx[{dx(j>)'^ /2 + V{(f>)]. In higher dimensions, 
this equation can be reduced to the standard model of 
curvature-driven growth , whereby the domain walls 
move with a velocity proportional to their local curva- 
ture. In one dimension, of course, the domain walls are 
points, and domain coarsening is driven by the exponen- 
tially weak forces between adjacent walls, mediated by 
the exponential tails of the domain-wall profiles. 

For the case of interest, where the typical domain size 
L is large compared to the width ^ of a wall, the closest 
pair of domain walls annihilate while the remaining walls 
hardly move at all, leading to the following simplified 
model 1^ . The smallest domain is combined with its two 
neighbors to form a new single domain. This process is 
then repeated. Eventually the system reaches a scaling 
state in which the distribution P[l,a) of domain sizes I, 
when the smallest domain size is a, has the scaling form 
P{l,a) = a-^f{l/a). 

An important feature of this process is that if the do- 
main sizes are initially uncorrelated they remain uncor- 
related, because the merging of three domains into a sin- 
gle domain introduces no correlations HJl^. It is then 
straightforward to derive an equation of motion for the 
evolution of the distribution P{l,a) as a increases [p[-p^: 

dP/da = -P{a, a)5{l - a) + e{l - 3a) P{a, a) x 

nl-2a 

X / P{l',a)P{l-l' -a,a)dl'. (24) 

J a 

Introducing the Laplace transform with respect to 
(j){p,a) — Jg P(Z, a) exp(— p/) d/, gives 

d<f)/da = -P{a, a) exp(-pa) (1 - <f). (25) 

This equation can, in fact, be integrated exactly but, in 
the spirit of the previous sections, we shall first look for 
solutions of the form 

(f){p, a) = 00 (pa) + <f>i{pa) + ■■■ (26) 
P(a,a) = i/o(l) + -^/i(l) + ---, (27) 
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where fo{x), fi{x) are the scahng and correction-to- 
scahng functions for the domain-size distribution, (f)i 
are the corresponding Laplace transforms, and u is the 
correction-to-scaling exponent as usual. Inserting (|2^ ) 
and (|2^) into (p5|), and equating leading and subleading 
terms gives 



— = |- + 2/o(l)^^0o|0i 

. ,,, exp(-5) 2^ 

- /i(i) — : — (1 - n), 



(28) 



(29) 



where s = pa. With the boundary conditions (/>o(oo) = 
(which follows from the definition of 0o) and '/>i(0) = 
(which follows from the normalization of P{l,a)), these 
equations can be integrated to give 



(f>o{s) = tanh 



Ml) J exp(-i) 



(/.i(s) = /i(l)5-(l-02(,)) 



dt 



(30) 



exp(-f). (31) 



Consider first the scaling function 4'o{s). For small 

s, one obtains 0o(s) = 1 - 2e^-^o(^)^ s^-^^f^^ H , where 

7 = 0.577 • • ■ is Euler's constant. But, from the defi- 
nition of 4>o{s) as the Laplace transform of fo{x), one 
also has (f>o{s) = 1 — s{x)q + • • •, where {x)o is the mean 
domain size (in units of the minimum domain size a) 
for the scaling distribution. Comparing these expansions 
gives /o(l) < 1/2, with /o(l) = 1/2 when the first mo- 
ment of fo{x) exists, and /o(l) < 1/2 when it doesn't. 
We shall primarily consider the former case, appropriate 
to an initial distribution with a finite first moment. 

Now consider the small-s behavior of the correction- 
to-scaling function (t>i{s). The factor (1 — 0o(*)) ™ (0) 
behaves as s^-^o^^) for s — > 0, i.e. as s for /o(l) — 1/2. For 
this case, the small-s expansion of (f>i{s) contains a non- 
analytic term 5^+"^ for any uj > (with a In s factor when 
uj is an integer). This in turn implies that the correction- 
to-scaling function fi{x) has the power-law tail x~^^^^\ 
This means that if the domain-size distribution is given 
by the scaling distribution foix) plus a small perturba- 
tion with a power-law tail x~^'^^'^^ , the perturbation is ir- 
relevant (i.e. decays to zero) for uj > 0, and the correction- 
to-scaling exponent is uj. More generally, any initial dis- 
tribution with a power-law tail x~^'^^^\ with w > 0, 
flows to the scaling distribution /o(x) with /o(l) = 1/2. 
For < 0, the perturbation is relevant. In this connec- 
tion one should note that the scaling distributions with 
/o(l) < 1/2 have the power-law tail a;-[i+2/o(i)] ^ i.e. they 
decay more slowly than 1/x^. So any initial condition 
that falls off more quickly than one of these 'special' fixed 
distributions flows to the 'generic' fixed distribution with 
/o(l) = 1/2 (which may easily be shown to have an ex- 
ponential tail). 

We have shown that any finite correction-to-scaling ex- 
ponent UJ corresponds to a correction-to-scaling function 



with a power-law tail. This implies that corrections to 
scaling that decay faster than any power-law in the scal- 
ing variable x have an infinite correction to scaling expo- 
nent, i.e. that the amplitude of the correction does not 
have the power-law form a^^ assumed in (p6|). To clarify 
this case, we construct the exact solution of (|2^) for an 
arbitrary initial condition: 



0(p, a) — tanh 



da' P{a' ,a') exp(— pa') 



(32) 



The arbitrary function of p that formally appears inside 
the argument of the tanh, as an additive integration con- 
stant, must vanish because (f>{p, oo) = 0. The function 
P{a, a) appearing in ( ^2| ) is determined by the initial 
condition, for which a = oq, through the inverse Laplace 
transform 



P(a,a) 



c 



dp 
Ani 



exp(pa) In 



1 + (t>iP, «o) 
1 - Hp, ao) 



(33) 



where the contour C runs parallel to the imaginary axis, 
to the right of any singularities of the integrand. 

Equations ( p^ and (^3|) represent an exact closed-form 
solution for the full evolution from an arbitrary initial 
condition. For explicit results, however, one needs to be 
able to evaluate the Bromwich integral (|3^). For illustra- 
tive purposes, we consider two simple initial conditions. 
The first is the simple exponential distribution P(/, 0) = 
pexp(—fil), with Laplace transform 0(^,0) = fi/{p + fi) 
(i.e. oq — here). Inserting this in (p3) gives 



P{a,a) ;^[1 - exp(-2/ia)]. 
Za 



(34) 



One sees that P{a,a) approaches the scaling limit l/2a 
exponentially fast as a function of a, i.e. faster than any 
power. This is in accord with our result that a power-law 
correction to scaling is associated with power-law tails in 
the initial domain-size distribution. Inserting (|3^ ) into 
(p3), and expanding for fia ^ 1, gives 



(j)ip,a) = (j)o{s) 



^ — 2fia — s 



{I -Ms? 



(35) 



where 00 (s) (with s = pa) is the scaling function. Again, 
the correction is exponentially small. 

Our second example shows that the correction to scal- 
ing can be an oscillatory function of a. For P{1,0) ~ 
/i^/ exp(— uT), one obtains (?!)(|5, 0) = /i^/(p+/i)^. Inserting 
this into (|3^ ) gives P(a, a) = [1 — 2 cos(^a) exp(— ^a) + 
exp(— 2/ia)]/2a, so that the leading correction to scaling 
is the term involving cos(/ia). The difference between this 
case and the first example is that in the former the singu- 
larities of the function h{p) = \n{[l-\-(j){p, O)]/[l-0(p, 0)]} 
appearing in ( |3^ ) are branch points on the real axis at 
p — and — 2/i, whereas in the latter they are branch 
points at p = 0, — (1 ± i)fi and —2^. The leading scaling 
part l/2a in P(a, a) is associated with the singularity of 
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h{p) at p = 0. The leading correction to scaling derives 
from the singularity with the largest real part (exclud- 
ing p — 0). If this singularity is off the real axis, the 
leading correction to scaling will be oscillatory. In all 
the other models considered in this paper, however, the 
leading correction to scaling is monotonically decreasing 
function of time. 



V. THE OJK THEORY 



7r 



d-i 



7o + (^^+4j 70+470 = (38) 

x\ I / d + Lu\ b , , ^ 

+ 4 j 7i + [-^) 71 + ^^I'o = . (39) 



In fact these equations are identical to equations ( |l9|) and 
( pO| ) of the large-n theory, with the same boundary con- 
ditions, and the results for 70 and 71 can be simply read 
off: 



The theory of Ohta, Jasnow, and Kawasaki (OJK) g] 
for the pair correlation function of a nonconserved scalar 
field starts from the Allen-Cahn equation relating the in- 
terface velocity to its local curvature ||^. The theory is 
expressed in terms of a 'smooth' auxiliary field m(x, t) 
whose zeros give the positions of the interfaces between 
domains of the two phases. The field m obeys the dif- 
fusion equation dm/dt — V^m (we absorb the diffusion 
constant into the timescale). The initial conditions are 
taken to be is taken to be gaussian. The normalized pair 
correlation function of m is 



7(r, t) 



(to(x, t) to(x + r, t)) 
(to(x, t)2) 



(36) 



while the correlation function for the order-parameter 
field 6 is 



C(r, t) = (sgn[rn(x, t)] sgn[m(x -|- r, t)]) = — sin ^ 7(r, t) 

TT 



(37) 



Note that by using the sgn function to relate the order 
parameter to the auxiliary field we are working in the 
'thin wall' limit. This means that we are neglecting any 
corrections to scaling associated with the finite width of 
the domain walls, and are focusing instead on corrections 
deriving from the initial conditions. This is the same ap- 
proach that we are using throughout this paper. 

It is simplest, in the first instance, to compute the cor- 
rections to scaling associated with the function 7(r, t). If 
we define h{r,t) — (to(x, t) m(x + r, i)), then '^{r,t) — 
h{r,t)/h{0,t). From the diffusion equation for m, one 
obtains immediately {l/2)dth = V^h, giving (1/2)9(7 — 
V^-f + a{t)-/, where a{t) = ~[dt\nh{0,t)]/2. Since 
h satisfies a diffusion equation, we know that for a 
short-ranged initial condition, peaked around r — 0, 
h{0, t) '-^ t^'^/'^ asymptotically, giving a{t) d/it for 
large t. It is convenient, however, to incorporate cor- 
rections to scaling in L{t) through the requirement that 
a{t) = A/2L2 exactly, i.e. (1/2)9*7 = + i\/2L^)-/. 
Choosing the scale of L such that L t^^^ asymptoti- 
cally then fixes A = d/2, as in the large-n theory. Writing 

7(r,t) = 7o(r/i)+i-'^7i(r/L), audi = l/2L-h&/i^+'^, 
as usual, substituting into the equation for 7, and equat- 
ing leading and next-to- leading powers of L, gives 



71 (a;) 



7o(x) = exp(-a;V8), 
[bx'^/32{d+2)] exp{-x^/8). 



(40) 
(41) 



The correction-to-scaling exponent is w = 4, as in the 
large-n theory. 

The order-parameter correlation function C (r, t) has 
the expansion C(r, t) ^fo{r/L)+L~'^fi{r/L) + ---. The 
functions fo{x) and fi{x) follow immediately from (|37 

m and my. 



fo{x) ^ - sin ^70 (a;), 

TT 



nix) 



(42) 
(43) 



Note that /i(x) ~ x^ for small x, so that the linear term 
in x from /o is not modified - our definition of L has 
ensured that areal density of domain walls, and therefore 
the energy, varies as 1/L exactly. 



VI. THE MAZENKO THEORY 



An alternative 'gaussian closure' theory to that of OJK 
has been proposed by Mazenko For nonconserved 
scalar fields, the pair correlation function C{r, t) satisfies 
the closed equation 



\dC_ 
2~dt 



1 



nSoit) 



tan 



(I 



C 



(44) 



The function S'o(t) is defined as (m(r,t)^), where m is a 
gaussian auxiliary field as before. For present purposes, 
however, it is sufficient to note that Sq has dimensions 
(length)^. In fact, it is convenient to define the coarsen- 
ing length scale L{t) by Sq = L'^/X, where A is a constant 
whose value is fixed by physical requirements |3|. This 
definition of L is in accord with our previous definitions, 
as we shall see. 
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FIG. 1. The scaling function fo{x) in the Mazenko theory 
for d = 2. 



FIG. 2. The correction-to-scaling function fi{x) in the 
Mazenko theory for d — 2. 



Writing 5*0 = L^/\ in (||), setting C{r,t) = fn{r/L) + 
L-^fi{r/L) + ■ • •, dL/dt = 1/2L + b/L^+'^ + and 
equating leading and next-to-leading powers of L in the 
usual way gives the following equations for the functions 
fo{x) and /i(x): 

+|/. + 5i/; = 0. (46) 



Eq. (^ ) is to be solved with boundary conditions /o (0) = 
1, /'(O) = -(l/7r)[2A/(d-l)]i/2, the latter following from 
a small- a; analysis of the equation. The parameter A is 
determined as follows For large x, ( |45| ) reduces to the 
linear equation + [{d -l)/x + + (A/2)/o = 0. 

This equation has the same form as the large-n equa- 
tion (|l^), and the two linearly independent large- 2; so- 
lutions have the large- a; forms x~'^^ and ~ exp(— x^/S) 
discussed in section III. If the equation is integrated for- 
ward from X — 0, for general A a linear combination of 
the two large-x solutions will be obtained. But a power- 
law decay is unphysical if the initial conditions contain 
only short-range spatial correlations. The parameter A is 
fixed by the requirement that this unphysical power-law 
decay be absent. Imposing this condition through a nu- 
merical solution of ( [l5| ) gives A — 0.711277 for d ~ 2, and 
A = 1.327411 for d — 3 (correct to the number of figures 
given). The exponent A, which describes the decay of the 
autocorrelation function via {(j){x,t)ip{x,0)) ^ L{t)~^, is 
given by A = d — A within this theory For d ^ 00, 
the OJK solution is recovered and A — > d/2. 



The correction-to-scaling uj is determined in a similar 
way from Eq. (HQ). The boundary conditions on this 
equation are /i(0) = = /{(O). For large-x, the ho- 
mogeneous part of this equation (i.e. without the final 
term) reduces to the homogeneous part of the equiva- 
lent large-n equation (po|). The two linearly independent 
large- a; solutions decay as ^ exp(— a;^/8) and ~ a;^^'^^". 
The exponent omega is determined by the condition that 
the power-law decaying term is absent from the large- a: 
solution. Implementing this condition numerically gives 
Lu = 3.8836 for d = 2 and cu = 3.9030 for d = 3. For 
d 00, the OJK result w = 4 is recovered. 

Figures 1 and 2 show the scaling and correction-to- 
scaling functions for d = 2 (those for d = 3 are very 
similar). The amplitude of the correction to scaling func- 
tion is fixed by the constant b in equation ( |4^ ) . The value 
b = 2 was used in Figure 2. 

Note that, within the Mazenko approach, the expo- 
nent u is nontrivial. Indeed in this theory to is on the 
same footing as A, which is known to be nontrivial in 
general. This suggests that, like A, the correction-to- 
scaling exponent uj is in general a nontrivial exponent of 
phase-ordering kinetics. 



VII. CONCLUSION 

An understanding of the corrections to scaling is im- 
portant in analysing data from experiments or simula- 
tions. In this paper we have computed the correction-to- 
scaling exponents to, and the corresponding correction- 
to-scaling functions, for a number of models of phase- 
ordering dynamics. For simplicity we have focussed on 
the correction to scaling associated with the approach 
of the order-parameter morphology to its scaling form, 
i.e. the correction arising from a nonscaling initial condi- 
tion, and have suppressed contributions associated with, 
for example, the finite thickness of the domain walls by 
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working in the thin- wall limit. 

In simple soluble models such as the ID Glauber model, 
the n = 00 model, and the OJK theory, the exponent w 
takes an integer value, a> = 4. Within the Mazenko the- 
ory, however, oj takes nontrivial values, suggesting that 
this exponent is, in general, a new nontrivial exponent 
of phase-ordering kinetics. The ID TDGL equation is 
anomalous in that corrections to scaling vanish exponen- 
tially fast as a function of the smallest domain size. 

It is important to note that the values of w computed 
in this paper correspond to a special, though natural, 
definition of the domain scale L through the mean en- 
ergy density. Specifically the energy density is defined to 
be exactly 1/L (up to an overall constant) for the scalar 
models, i.e. the areal density of domain walls is propor- 
tional to 1/L, except for the ID TDGL equation where 
the scale length was conveniently chosen to be the size 
of the smallest domain. For the large-n non-linear sigma 
model, the energy density is defined to be exactly l/L"^, 
again up to an overall constant. These definitions of L are 
optimal in that they minimize the corrections to scaling, 
expressed in terms of L, i.e. they maximize oj. 

In the present paper wo have considered systems with 
nonconserved order parameter only. In future work we 



hope to discuss corrections to scaling in systems with 
conserved order parameter. 
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